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Abstract 

We present the solution to London's equations for the magnetic fields of 
a vortex oriented parallel to the planes, and normal to a crystal face, of a 
layered superconductor. These expressions account for flux spreading at the 
superconducting surface, which can change the apparent size of the vortex 
along the planes by as much as 30%. We compare these expressions with 
experimental results. 



A. Introduction 

Recently, scanning SQUID microscope magnetic imaging of interlayer vortices trapped 
between the planes of layered superconductors has been used to make direct measurements 
of the interlayer penetration depth in several layered superconductors |], @, @], [0- These 
experiments provide local measurements of the interlayer supercurrent density, which have 



implications for the validity of the interlayer tunneling model as a candidate mechanism 
for superconductivity in the high critical temperature cuprate superconductors. 

To date the quantitative modelling of these experiments has assumed that the vortex 
fields at the superconductor-vacuum interface are the same as those in the bulk, neglecting 
the well known effect that the magnetic fields from vortices spread as they approach the 
superconductor-vacuum surface from within the superconductor. Exact theoretical expres- 
sions exist for a vortex in an isotropic London's model ||. For a vortex oriented perpen- 
dicular to the surface in a superconductor with an isotropic penetration depth A, the fields 
above the surface can be approximated by a magnetic monopole located a distance A below 
the surface ||, [J7|. This means that the spatial extent of the magnetic fields at the surface 
is larger than in the bulk of the superconductor. If the bulk expressions were used to fit 
data at the surface, the fitted value of the penetration depth would be longer than the real 
value. This effect must be accounted for in making quantitative estimates of the penetration 
depths by magnetic imaging measurements. For this purpose it is useful to examine vortex 
spreading at the surface for a highly anisotropic superconductor, since recent experiments 
have studied vortices in superconductors with \ c /\ a b ~ 10-100. 

It is well known that the anisotropic London model is appropriate for describing a stack 
of Josephson-coupled superconducting layers at length scales large compared to the in- 
terlayer spacing. In this paper we present an exact solution of London's equations for a 
straight vortex approaching a superconductor- vacuum interface normal to the interface, in 
an anisotropic superconductor. We show how flux spreading near this interface effects the 
magnetic fields above the interface, and show that there is good agreement between these 
theoretical results and scanning SQUID microscope measurements on single crystals of the 
layered high-T c cuprate superconductor Tl 2 Ba2Cu0 6+( 5 (Tl-2201) 
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B. The model 



A method for finding the field distribution of a straight vortex crossing a plane surface of 
an anisotropic superconductor has been developed in Ref. H. We will outline this method 
and apply it to the case of a vortex, oriented along b in the ab plane of a uniaxial material, 
which crosses the plane face ca of the crystal. For a vortex not too close to the crystal 
corners, the crystal surface ca can be taken as an infinite plane. We choose the coordinates 
x,y,z corresponding to c,a,b of the crystal as shown in Fig. 1. Then the mass tensor is 
diagonal: m xx = m 3 ,m y y = m zz = rri\. The standard normalization m^m 3 = 1 is implied. 
The method consists of solving London's equations for the field inside the superconductor 
and matching the result to a solution of Maxwell's equations in the vacuum outside the 
sample. For the isotropic case, the problem is simplified by the cylindrical symmetry of the 
field distribution ||||. This is not the case for anisotropic materials, and a more general 
approach is needed. 

Inside the superconductor, the field h(r, z), with r = {x,y}, satisfies the London equa- 
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hi - —\ 2 m k ieu s ^p- = (f) 5(r)5 iz . (1) 
c ox s 

Here, j is the current density, O — hc/2e is the superconducting flux quantum, and the 
average penetration depth A = (A^Ac) 1 / 3 (\ 2 ab = mi A 2 , A 2 , = m^\ 2 ). 

Deep inside the superconductor, the field h(r) has only a z component. However, near 
the exit from the sample at z — 0, the vortex "opens up" and h x , h y are no longer zero. In 
other words, Eq. ([!]) is a system of three linear differential equations for h x ,h y , and h z with 
a non-zero right-hand side (RHS). The general solution is then 

h = h (0) + h w , (2) 

where h^ ) solves the homogeneous system with zero RHS, whereas hw is a particular 
solution of the full system ([I]). The latter can be taken as the field of an infinitely long 
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unperturbed vortex along z; this assures correct singular behavior at the vortex axis. The 
Fourier transform of this field is 

With this choice of , the field h^ ) is the correction due to the surface of the unperturbed 
vortex field h^. We note that the Clem-Coffey result for a vortex parallel to the layers 
of a Josephson coupled layered superconductor reduces to if one disregards the core 
correction [ |TTf . 

Because the only sample boundary is parallel to the plane xy, we Fourier transform 
Eq. (|ip with respect to x, y. We are then left with the system of equations for h(k, z) = 
J dr exp(— i\t ■ r)h(r, z): 

mih" — (1 + miky)h x + m\k x k y h y — imik x h' z = , 

mik x kyh x + m 3 h'' — (1 + mxk^.)h y — im 3 k y h' z = , (4) 
imik x h' x + im 3 k y h' y + (1 + m\k\ + msky)h z = (f) . 

For brevity, we have set the average A as the unit of length so that \ 2 ab = A 2 m x and \ 2 C = \ 2 m 3 
are replaced with mi and m 3 ; the prime in the above equations denotes d/dz. The field 
(k, z) satisfies the homogeneous system of linear second order ordinary differential (with 
respect to the variable z) equations, i.e., it is a linear combination of exponential functions 
of z: 

h (0) (k,^) = ^H (n W. (5) 

n 

The z independent coefficients H^(k) and a n (k) are still to be determined. Each term in 
the sum (|^) should satisfy separately the system @ with zero RHS. Omitting the label n 
we write this system as 



A tJ H 3 = (6) 



with a symmetric matrix A«: 



A xx = 1 + mik y - mi a , 

A xz = im 1 k x a, (7) 

A TO = 1 + mi^ - m 3 a 2 , 

A y2 = im 3 k y a , 

A zz = 1 + mi^ + m 3 /Cy . 

The determinant of this matrix must be zero, which provides all possible values of a: 

/l + miAr\i/2 

"1,2 = ±( ) , (8) 

K mi J 

A + mikl + m 3 k 2 \ 1/2 
«3,4 = ± • (9) 



m 3 

Deep inside the superconductor, the surface correction hS°^(z — > —00) must vanish, implying 
aii and a 3 must be positive. The homogeneous system (|) allows one to express (for each 
of these a's) two out of three components Hi in terms of the third. We obtain after simple 
algebra: 

ffM = ii±^ff«, if«=z^tf«; (10) 
m\k x a\ ol\ 

H x 3) = 0, H®=if.H,(3). (11) 

Ky 

Thus, the field inside the sample will be determined completely after Hf' and are found 
from the boundary conditions at the sample surface. 

The field outside the sample is described by divh = and curlh = 0, so that one looks 
for h = V93 with V 2 (p = 0. The general solution of Laplace's equation which vanishes at 
z — > 00 is 

The 2D Fourier transform is defined by 

<p(k) = e kz J d 2 r ip(r, z) e" ik r . (13) 



The boundary conditions at the free surface z = consist of continuity of the three field 
components: 

ik x <p = + #00 , 

ik yV = H^ + Hf\ (14) 
-ky = h^+H^ + H {3) . 

The components H^ y 3 ^ are expressed in terms of H^ 1 ' 3 ^ in Eqs. (]10f) , (0), so that the system 
(14]) can be solved to find ip along with all H^s. We are interested here primarily in the field 
outside the sample: 

m 3 a 3 [mj^ a 3 (fc + a x ) + /ca 3 + fcjj] ' 

It is readily verified that for the isotropic material Eq. flT5|) reduces to the known result by 
Pearl: ip = —(f> /ai S k(k + a.i s ) where oti S is the isotropic version of either a% or 03 0. 
Since outside the sample h = V<£>, we have 

h x , y (k) = ik x , y (p(k) , h z (k) = -k(p(k) . (16) 

Then, for example, the 

h z (r,z) = -J-^ 2 k V (k)e^- k *. (17) 

In particular, the total magnetic flux through any plane z = zq is given by h z (k = 0, zq) = <po 
as expected. 

The field inside the sample is given by Eqs. (fj), (§). The coefficients H^ 1 ' 3 ) in 

h (0) (k, z) = H (1) e ai2 + H (3) e° 32 (18) 

are obtained by solving (|I4]) , (|TD|), and (PI): 

H«') = V (k)[ ife ,,-^%,-^%}, (19) 
L l+m-\k£ l + m^ki> 



+ niikl l + mik% 
iky ky 
+ mi^ ' a 3 (l + miA;2) 



H,3) ^( k H°.T^-^rii!)}- < 20 > 



Thus, for example, 



hJr, z) 



d 2 k 

WT 2 



ik x ip(k) e 



(21) 



For what follows, we will only concern ourselves with the fields outside of the supercon- 
ductor. In the experiment of Ref. |l|, the component h z was probed with a SQUID pickup 
loop which was much larger than the penetration depth A a ft. One therefore expects the 
instrument to measure a flux nearly equal to the pickup loop size times 



T-Cz{x,y) 



h z (x, y, z) dx . 



(22) 



The vortex spreads as it approaches from below the superconducting surface in the re- 
direction as well as in the y-direction; nevertheless numerical estimates show that under 
typical conditions the experimental signal is well represented by Eq. (p2[). Then we obtain: 



elk-). 



2ir m 3 a (a + \k y \) 



(23) 



with a = a^(k x = 0) = (m 3 1 + k 2 ) 1 ^ 2 . In conventional units, 



7-Lz(y,z) 



l-oo 2tt a (a + X c \k y \) 

with a = (1 + \ 2 k 2 ) 1 / 2 . It is worth noting that the quantity H z (y, z) depends only on A c . 
After the substitution X c k y = sinhu, Eq. (p3|) takes the form 
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7rA, 



n z (x,y) 



du e 



Re 
Re 



o 
rl 

w 



-u-z'sinhu cos ( y ' sinhu ) 
due -u- WS inhu 

7r/ E 1 (w)+Y 1 (w) 



(25) 



Here, y' = y/X c , z' = z/X c , and w = (z + iy)/X c ; Ei and Y\ are Weber's and Neumann 
functions, see Ref. HI! 
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For \w\ <C 1 (both y and z are small relative to A c ) we have [12 
H z (y,z 



Re 



z 



W /■ w 1 
1 + — I In — + 7 

2 V 2 ' 2 



(26) 



. Jy 2 + z 2 



(7 = 0.577 is Euler's constant). If \w\ 3> 1 (at least one of y or z is large relative to A c ), we 
obtain: 



W.(V>*)_oJ-l _i_ +0(H -4 



7rA,. — z — - — = Re 



^A c (^ 2 -i/)A 2 



+ .... (27) 



I/ 2 + z 2 (z 2 + y 2 ) 2 

At large distances from the vortex exit, the second term can be neglected and A c drops out 
of the result; this is expected since the field there is approaching the Coulomb form with no 
trace of material properties. 

To form a complete picture of the field distribution outside the sample, we imagine that 
the same SQUID probe is oriented in the xz plane so that the y-component of the field 
(integrated over the probe area) is measured. Then the SQUID flux will be nearly equal to 
the pickup loop size times 



H y (y,z) = / h y (x,y,z)dx 

J — oo 

/oo A U 
^ik y <p(0, k y )e ik y y - k y z . (28) 
-oo 2n 

The two-dimensional field 7i satisfies divli. = curYH = 0. Hence, it can be written as 
7i = V$ with V 2 $ = 0. This implies that 7i z and 7i y are real and imaginary parts of the 
same analytic function given in Eq. (p5[) : 

H y (y, z) = ^Im{- + \ feH + Y x {w)\ } . (29) 



7rA c 2 

In particular, the asymptotic form of this function for w ^> 1 can be written as 

1 1 1 



(30) 



w w 2 w + 1 ' 

which implies that at large distances the field TC behaves as a field of a 2D "charge" situated 
at z — — A c . 



C. Results 



Figure 2(a) shows streamlines of 7i(y, z) for a single anisotropic vortex centered at x — 0, 
y = 0. These streamlines were generated numerically as follows: the starting points of the 
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lines were at z/\ c =-3, with a spacing in y between the lines proportional to (dHz/dy)^ 1 
at z/X c = —3. Small steps z = z + Az, y = y + Ay, with Az = 5 sin 0, Ay = 5 cos 9, 
were generated, with # = tan _1 (7i 2 (?/, z)/TC y (y, z)). The fields were recalculated at the new 
positions, and then the process was repeated until | z/X c |> 3 or | y/X c |> 3. Figure 2(b) 
shows the results if field spreading below the surface is neglected: 

b z {x, y,z<0) = -^-K (R), (31) 

where K is a modified Bessel function of the second kind of order 0, R = [(s/2X a b) 2 + 
(x/X a b) 2 + (y/Ac) 2 ] 1 / 2 , and s is the interplanar spacing [|TTJ. For s <C X ab , Eq. ( |31| ) has the 



Fourier transform given in Eq. (Q). The fields for z > are given by 

d 2 k 



b z (r,z) = J j^~M k ) ^' r ~ kZ » ( 32 ) 
where 

b z (k) = J d 2 kb z (x, y, z = 0)e- lk r . (33) 

TM. The fields in the x and y directions are treated similarly, using the relations b x (k) = 
—ik x b z (k)/k, and ^(k) = —ik y b z (k)/k. The calculated fields were integrated over x and 
streamlines were generated just as for the exact London expressions. This is the procedure 
used to model the experimental results in Ref. In both cases the fields extend into the 
vacuum nearly isotropically at large distances, as if a point monopole source were placed 
near z = —X c . In the full treatment, Figure 2(a), the fields spread as the vortex approaches 
the surface from inside the superconductor. 

A one- dimensional rendering of our results above the surface of the superconductor is 
shown in Figure 3, which plots TrX c H z (y, z)/4> as a function of y/X c for several values of 
z/X c . For comparison, the results neglecting vortex spreading are also shown. Note that 
7i 2 (0,0), which approximately indicates the peak signal in an experiment, is overestimated 
by a factor of tt/2 if vortex spreading is neglected. Also, the full width at half-maximum of 
the flux contour is 1.87A C for the full theory, while it is 1.37A C if flux spreading is neglected. 
Therefore the neglect of flux spreading could result in an overestimate of A c by 30%. 
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D. Comparison with experiment 



Previous analyses of experimental data 0,0, neglected the effect of vortex spreading at 
the surface. In retrospect, this neglect was not unreasonable, given the quantitative agree- 
ment between the interlayer coupling strength obtained from vortex imaging measurements 
and from the Josephson plasma resonance 0. With the full theory presented in this paper, 
it is now possible to quantitatively examine this assumption. However, we note that there 
are additional systematic experimental uncertainties in this technique. These errors include 
the effect of macroscopic screening currents which may create a slightly inhomogeneous 
background, a relative angle of 10-20 degrees between the surface of the superconductor and 
the SQUID pickup loop, the uncertainty in the exact value of the height of the pickup loop, 
and some effect of the leads to the pickup loop on the effective shape of the pickup loop. 
We estimate that these errors may also be as large as 30%. 

Figure 4 shows grey-scale images of 5 interlayer vortices in Tl-2201 magnetically imaged 
using a scanning SQUID microscope ||]. In these experiments the experimental signal is 
equal to the integrated magnetic flux through the pickup loop. The vortices appear elliptical 
in shape, with the long axis (parallel to the planes) nearly vertical in these images. The 
spatial extent of the vortex images perpendicular to the planes is limited simply by the 
size of the pickup loop. The spatial extent parallel to the planes is set primarily by the 
interplanar penetration depth. Figure 5 shows cross-sections through the experimental data 
along the direction parallel to the planes as indicated by the dashed lines in Figure 4. 

To generate a theoretical expression for fitting the experiment results, we use the full 



expression (Eq. (lq)) for the magnetic fields, taking the a-axis penetration depth equal to 
0.17/mi . An evaluation of Eq. ([15]) gives the z-component of the field at a given height 
Zo above the sample surface. This field is then numerically integrated over the shape of 
the pickup loop to obtain the total theoretical flux <& s (x, y, z) through the pickup loop as a 
function of the pickup loop position and the c-axis penetration depth A c . In this case the 
pickup loop used was a square 8.2/mi on a side, with a 1/xm linewidth, and a superconducting 
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shield 5/xm wide which extends to the top corner of the pickup loop, as indicated by the 
inset of Figure 4. For our modelling we add to the flux through the pickup loop one third 
of the flux intercepting an area 5/im by 5/im on a side, starting at the upper corner of the 
pickup loop, to account for flux focussing effects from the superconducting shield. The solid 
lines in Figure 5 are best fits of the cross section § s (x — 0,y, z — z ) to the experimental 
data, using the interlayer penetration A c , the height of the pickup loop z , and an offset flux, 
representing a small uncertainty of the background signal, as free parameters. The best fit 
values for each vortex are displayed in the figure, with uncertainties assigned using a doubling 
of the x 2 value as a criterion. Reasonable agreement is obtained between experimental and 
theoretical cross-sections. The average value for A c , using a weighting inversely proportional 
to the square of the uncertainties, is A c = 18.3 ± 3/im. 

For comparison, a similar analysis of the same vortices neglecting the effect of vortex 
spreading resulted in a weighted average best-fit value of A c = 19 ± 2/im. This comparison 
is surprising at first, since the theoretical FWHM at the surface is reduced by 30% when the 
vortex spreading is neglected. However, the uncertain value of the height of the pickup loop 
compensates for the assumption of negligible vortex spreading below the surface. When the 
spreading below the surface is neglected, the fitting routine compensates by picking a higher 
value of zq, thereby moving the spreading to the vacuum rather than the superconductor. 

In conclusion, we have presented a solution to London's equation for the case of an inter- 
layer vortex approaching a superconducting surface normal to the surface and parallel to the 
planes. This model is appropriate for experiments which magnetically image interlayer vor- 
tices. Good agreement with available experiments is obtained with this model, allowing the 
quantitative determination of the interlayer penetration depths from these measurements. 

We would like to thank D.G. Hinks, T.W. Li, and Ming Xu for supplying the Tl-2201 
crystals used for the SQUID images shown in this paper. We would also like to thank M.B. 
Ketchen for the design, and M. Bhushan for the fabrication, of the SQUIDs used here. 
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Figure Captions 



1. Geometry and axes used in this paper. A single vortex, centered at x = 0,y = 0, 
emerges normally to the ac face of the crystal located at z — 0. 

2. Stream-line mapping of H(y,z), the magnetic field integrated over x, for a single 
inter layer vortex shown in Fig. 1. The spacing of the stream lines is proportional to 
(dHz/dy)^ 1 at z — — 3A C . Fig. 2a is the present model, which includes vortex spreading; 
Fig. 2b is for a model which neglects field spreading below the surface. 

3. Plot of H z (y,z) for a single interlayer vortex normal to the superconducting surface 
z = 0, as a function of y/X c , for fixed z values, where y is the distance from the center 
of the vortex along a plane direction, and A c is the interlayer penetration depth. 

4. Grey-scale images of 5 interlayer vortices in the high-T c cuprate superconductor 
Tl-2201. The scaling corresponds to O.110 O (I), O.160 O (II), O.130 o (ni), O.220 O (IV), and 
O.190o(V) full-scale variation from black to white in the integrated flux through the 
pickup loop. 

5. Cross-sections along the plane directions through the images of Figure 4. Each suc- 
cesive curve is offset by 0.1 unit for clarity. The dots are the data; the lines are fits to 
the present model as described in the text. 



13 




-60 -40 -20 20 40 60 

Position(|im) 




20/xm 



